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Abstract 



The total and differential cross sections for the production of ^scb baryons in electron-positron 



^ ' collisions are calculated at the Z-boson pole. 

> 
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Oi 1- Introduction 

Q ■ Baryons involving two or three heavy quarks (c, b) have not yet been observed experimen- 

i-C ' tally. Theoretical investigations into the mass spectra of hadrons that contain two or more 

' ! heavy quarks, the cross sections for their production in various processes, and their lifetimes 
Mh. 

^ ■ and decay modes form a rather new line of research in particle physics. For an overview of 
these investigations, the interested reader is referred to [1]. Calculations presented in the 
^ . literature that are aimed at determining the production cross sections for baryons featuring 
d • two heavy quarks either rely on the fragmentation approach [2], or treat the production 

of unbound four-quark states in a bounded phase space [3], or consider (as is done in the 
majority of the most recent studies on this subject) the production of the relevant heavy 
diquark [4-9]. In particular, the production of doubly heavy baryons in electron-positron 
collisions was analyzed in [3, 9]. No attention has so far been given to the production of 
triply heavy baryons. 

Investigation of the mechanisms responsible for the production of multiply heavy hadrons 

is of interest from the theoretical point of view since this provides the possibility of further 

verifying QCD (more precisely, of obtaining deeper insight into this theory). This involves 

testing both its perturbative aspect, which is used to describe the simultaneous production 
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of a few quark pairs, and QCD-inspired non-perturbative models of bound states. We recall 
that, even in apparently obvious cases, the results of QCD calculations often unexpectedly 
prove to be at odds with experimental data, as was, for example, in the hadronic production 
of J/V' particles. On the other hand, a derivation of theoretical cross-section estimates 
is of importance from the point of view of applications to searches for such particles and 
investigation into their properties. The presence of two or more heavy quarks in a hadron 
substantially affects the properties of its weak decays. At the same time, the reliability of 
theoretical predictions is higher for such hadrons, and this makes it possible to test model 
concepts more thoroughly. 

In the present study, we reveal some features of Qsc& baryon production in electron- 
positron annihilation. The choice of a process where the initial state is purely leptonic was 
motivated, on one hand, by the fact that the relevant calculations are simpler here than in 
the case of hadronic production and, on the other hand, by the fact that the case of leptonic 
production offers a number of advantages in what is concerned with a possible experimental 
observation, which include favorable background conditions and a precisely known initial 
energy. 

At the quark level, the subprocess e^e+ sscchh which is of order a^a^ in conventional 
perturbation theory, is associated with the process being considered. In evaluating the square 
of the relevant matrix clement, use is made of the method that was proposed in [10] and 
which is referred to as the method of orthogonal amplitudes (previously, this method was 
employed, for example, in [7, 8]). The fusion of the product s, c, and h quarks into a Vtscb 
baryon is described within the standard nonrelativistic approximation [11-13]. A detailed 
account of the technical facet of our calculations is given in Section 2. The numerical results 
obtained on the basis of these calculations are discussed in Section 3. 



2. Computational procedure 

Our calculations are based on considering the partonic process 

e-(Pe-) + e+(Pe+) ^ + C(p2, C) + KP3, X) + S(P4, O + c(p5, C') + KP6, x'), (1) 

where the parentheses referring to the colliding particles enclose their 4-momenta, while the 
parentheses referring to the product quarks and antiquarks enclose their 4-momenta and color 
indices. As usual, we disregard Feynman diagrams featuring the electroweak interaction of 
the quarks involved. In other words, we consider only those diagrams where the interaction 
between the quarks is mediated by gluons. For this class of diagrams. Fig. 1 shows nine 
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basic diagrams in which quark and gluon Unes are connected in different ways. Here, six 
(three) nonequivalent dispositions of the ss, cc, and bb hnes correspond to each of the first 
seven (last two) diagrams in Fig. 1 - that is, diagrams 1-7 (diagrams 8 and 9). Considering 
that the annihilation channel may involve either a photon or a Z boson, we conclude that 
the total number of diagrams in question is 96. 

The matrix element for the process in (1) can be represented in the form 

" 24 cos^^ (« - M| + iMzTz) 6s ^ ^ ^€«'CC'xx' ' 

where 

aZ Tjssccbb 1 TDSsbbcc , jDCCssbb , jDCcbbss , r^bbsscc , r^bbccss /q"\ 

^S«'CC'xx' - -^^'CC'xx' + ^««'xx'CC' + ^CCa'xx' + ^CCxx'^' + ^xx'^C'CC + ^xx'CC'K'' 
^ll'cfxx' = {[(P2 +Pi +P4)' - m2]-i[(p,+ +p,- -p^f - ml]-\p^ +P4)-'x (4) 

X (pi + P2 + P4 + P5)~^u{P2)T^lTi(;'Y{P2 + Pl + pA + mc)-f5v{-p5) X 

X u{p3)T^^rj\pe+ + Pe- ' Pg + mb)je{9v " 9\l5)v{-pe) + 

+ [{P5+Pl +P4)' - m^]-'[(pe+ +Pe- - Psf - mD'^pi +P4)"' X 

X (pi +P2+P4+ P5)~^u{P2)T^lTt^>-fs{-p5 -P1-P4 + mc)YH-P5) X 

X u{p3)T^^,-fs{-pe+ - Pe- + p3 + mbh^gy - 9Al5)v{-pe) + 

+ [{P5+Pl+P4:f - ml]~'^[{pe+ +Pe- - P&f " ^1]'^ {pi +^4)"^ X 

X (pi + P2 + P4 + P5)"^'li(P2)7ciT;^'7<5(-P5 - Pl - P4 + mc)Yv{-p^) X 

X u{p3)T^^,-f\p^+ +pe- - P6 + rribheigv - QAlsM-pe) + 

+ [(P2 +Pl +P4)^ - ml]-^[{pe+ +Pe- - Psf " + Pa)'^ X 

X (pi + P2 + P4 + p^y^u{p2)T^iTi^,'^'' {p2 + Pl + P4 + "^c)7<5^(-P5) X 
X u{p3)T^^,^e{-pe+ - Pe- + P3 + mb)'y\gv - gA'y5)v{-p6) + 

+ [(P2 +P5 +P3)^ - ml]-^[{pe+ +pe- -pef - ml]-\pi +P4)"^ X 

X (P2 + P5)~MP2)T^C'^5v{-p5)u{P3)TilTil,-f\p2 + Ps + P3 + ^fe) X 

X Y{Pe+ +Pe- -P6 + mb)^e{9v ' fi'A75)^^(-P6) + 

+ [{P2 + P5 + PsT - ml]-^[{pi + p4 + pef - ml]-\pi + p^Y^ x 

X (P2 + P5)-'m(P2)7Jc'^5^(-P5)^(P3)^x^^/x'7^(^2 + Ps + Pa + rUb) x 

X 7£(-Pi - P4 - P6 + mb)Y{9v - 9a15)v{-P6) + 

+ [{Pi+PA+PeY -ml]-^[{pe++Pe- - Psf - ml]-\pi + p4)-^ x 

X (P2 + P5y^u{p2)T^(;,-fsv{-p5)u{p3)T^iT['^,%{-pe+ - Pe- + Ps + ^b) X 
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X l\-Pi -P4-P6 + mb)Y{9v - ^a75)^^(-P6) - 

- (V2)[(pe+ +Pe- -Pef - ml]-\pi+pi)-''{p2 X 

X (pi +P2+P4+ P5)-'r'l(-P? -pi+pi+ pDq"'' + 

X M(P3)T^x'75(Pe+ + Pe- " P6 + mb)-f^{gy - ^A75)^^(-P6) " 

- (i/2)[(pe+ +Pe- -P3)' - +p4)-2(p2 +p5)-2 x 
X (pi + P2 + P4 + P5)-V"'1(-P? - p! + + P^)^'^'^ + 

+ (-2p^ - 2p^)^'^^ + {2p'i + 2p'i)g''']u{^2)Tl^nA-^^) x 

X M(P3)T^^,7£(-Pe+ - Pe- +P3 + "^6)7^(5'^ " 9\lb)v{-V6)] X 
X '«(Pl)r^e7^f^(-P4)^^(-Pe+)7''(yy - 9Alb)u{^e-), 

and where the expression for -^^^^'c^c'xx' obtained from the expression A.^^ic^c^i^-xi by setting 
in it = 1, = 0, gy = Qg, and g\ = 0, with Qg being the electric charge of the 
quark q {q = s, c, b) in units of the electron charge e. Since the number of noncquivalcnt 
Feynman diagrams belonging to type 8 or 9 (Fig. 1) and differing from one another only by a 
permutation of the ss, cc, and bb lines is equal to three and since the quantity A^^,^^,^^, from 
Eq. (3) involves six terms that arc obtained from one another by permuting the quantum 
numbers of ss, cc, and bb quark-antiquark pairs, we introduce, in the eighth and the ninth 
term in expression (4), an additional (in relation to the Feynman formulation) factor of 1/2, 
whereupon we arrive at the correct results for this quantity and for the matrix element in 
(2). 

Let us consider in more detail the color structure of the matrix element given by Eqs. 
(2)- (4). Since any baryon is a color-singlet object, the scb state that is produced in process 
(1) must be an SU (3)c singlet that is contained in the tensor product of three SU (3)c triplets. 
It follows that the scb state must be fully antisymmetric in the color indices of the quarks; 
to take this into account, it is necessary to introduce, in the amplitude of the process, the 
antisymmetric tensor s^'^'^'/'n/S, which is normalized to unity. Since the initial electron- 
positron state is also a color singlet, the state of three unbound s, c and b antiquarks 
accompanying the product baryon must also be a singlet. We note in passing that the 
tensor e^''^T^^,T^iTi^,T^^, is already fully antisymmetric in its indices (' and x' therefore, 
it not necessary is s^'^'^'/v^ ^bat the projection operator be exphcitly present there. The 
presence of this operator is technically useful, however, since this makes it possible to perform 
summation over color indices at the amplitude level-that is, prior to squaring the amplitude. 
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Further, it is straightforward to prove the identity^ 

S^'^X^eCx' fabdj.a^^j.b^^j.d^^ = 0, (5) 

from which it follows that the contributions of diagrams involving a three-gluon vertex [eighth 
and ninth term in Eq. (4)] vanish. At the same time, all of the remaining terms in Eqs. 
(2)- (4) have the same color structure. Summation over color indices yields 



In performing summation over fermion polarizations (with the aid of the REDUCE system 
[14] for analytic calculations), we employed the method of orthogonal amplitudes. Briefly, 
the essence of the method is as follows. Suppose that we have the quantity u{p')Ru{p"), 
where u{p') and u{p") are spinors that obey the Dirac equation, while R is an operator 
that is expressed in terms of the matrices and their contractions with 4-vectors. In general, 
this quantity then admits a linear decomposition in terms of four orthogonal amplitudes 
wi = u{p')u{p"), W2 = u{p')Ku{p"), W3 = u{p')Qu{p"), and W4 = u{p')KQu(p"). That 
two different amplitudes are orthogonal implies the vanishing of the quantity obtained by 
summing, over the polarizations of the two spinors, the product of one of these amplitudes 
and the complex conjugate of the other. This is so if the 4-vectors K^^ and are orthogonal 
to the 4-momenta and p''* and p"^ to each other - that is K^p'^ — 0, K^jp"'^ — 0, Q^p''^ — 0, 
QnP"^ — 0, and KnQ^ — 0; otherwise, the 4-vectors and are arbitrary. 

Let us apply the method of orthogonal amphtudes to the specific problem at hand. For 
this purpose, we introduce the quantities 

Wsl = U{pi)v{-P4), Ws2 = u{pi)KsV{-p4:), (7) 

Ws3 = u{pi)QsV{-p4), Ws4 = u{pi)ksQsV{-p4), 

Wcl = U{p2)v{-P5), Wc2 = u{p2)KcV{-p5), 
Wc3 = u{p2)QcV{-p5), WcA = u{p2)KcQcV{-p->) , 



^Of three indices a, h and d corresponding to nonzero values of the structure constant Z"*"* two are 
always the numbers of Gell-Mann matrices such that they undergo no changes upon transposition, while 
the remaining index is associated with a Gell-Mann matrix that changes sign upon transposition. By 
simultaneously replacing the primed indices by unprimed ones and transposing the matrices T°, and T'^, 
we arrive at 

AiXfri'i'x T"- T'^ — —ACXo-i'i'x'n-'a rpb rpd _ __€'C'x'=.£Cx7io t'^ 

£ £ - £ £ - £ £ 

whence we immediately obtain relation (5). 
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where 



Wbl = ^i(P3)^^(-P6), Wb2 = u{p3)Kbv{-pe), 
Wb3 = u{P3)QbV{-p6): Wb4 = u{p3)KbQbV{-p6): 
Wei = v{-Pe+)KeU{Pe-), W^^ = ^(-Pe+ )Qe«(Pe- ) , 

K^^=e^-p-p^,p^^a,,, = £'^"'"^pi.P4pi^.a, (8) 

the 4- vectors a^o-, Oco-, Ofco-) a-nd Ogo- being arbitrary. 

Our problem is then described in terms of 128 orthogonal amplitudes of the form 

Wijkl = WsiWcjWbkWeh iJ: A) = 1, 2, 3, 4, / = 1, 2. (9) 

In order to find the coefficients Cijki in the expansion of the matrix element (2) in the 
amplitudes specified by Eq. (9), 

4 2 

M= Yl Y^ijkiWijkh (10) 

i,j,k=\ 1=1 

we multiply both sides of this equality by the quantity w*,j,j^,i,, sum the result over the 
polarizations of all fermions, and make use of the orthogonality of the different amplitudes 
Wijki- Denoting by \wijki\'^ the quantity obtained by summing, over the polarizations of all 
fermions, the squared modulus of the amplitude Wijki, we have 

Cijkl = { ■MKkl}/\^iJkl\'^- (11) 

polar. 

For the squared modulus of the relevant matrix element, summation over the polarizations 
of product particles and averaging over the polarizations of colliding particles is performed 
by the formula 

1 4 2 

1-^1' = 4 E E|c«iwlVu«l'- (12) 

i,j,k=l 1=1 

We note that we did not include the quantities v{—pe+)u{pf.-) and ■u(— Pe+)i^eQe^(Pe-) 
in the list of basic orthogonal amplitudes in (7), since the corresponding coefficients in the 
expansion of the matrix element [Eqs. (2)-(4)] vanish (this is because the traces that arise 
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in performing summation over the polarizations of massless electrons and positrons involve 
an odd number of Dirac 7 matrices). 

The question of why it is profitable to employ the method of orthogonal amplitudes is in 
order here. Upon directly squaring the matrix element specified by Eqs. (2)- (4), we would 
obtain, with allowance for the equality in (5), 3570 terms, and an individual operation of 
summation over particle polarizations would correspond to each of these terms. Within the 
method of orthogonal amplitudes, we compose one REDUCE code for evaluating traces and 
tensor contractions that corresponds to 84 terms in the quantity Adwim, whereupon we 
apply text editors (for example, joe or gedit) to perform obvious substitutions in this code, 
thereby obtaining REDUCE codes for evaluating all 128 coefficients Cijki- We note that 
analytic expressions for 128 coefficients Cijki occupy 370 Mb. 



3. Numerical results 

In order to describe a bound state of heavy quarks, wc make use of the nonrelativistic 
approximation [11-13], according to which the relative velocities of the quarks in a heavy 
hadron are assumed to be low. In the case of iS-wavc states, these velocities can be set to 
zero. Accordingly, the velocities of all three quarks in the final state of process (1) arc taken 
to be identical, while the momenta of the quarks are assumed to be proportional to their 
masses; that is, 

Pi = (ms/M)p, p2^{mc/M)p, p3^{mb/M)p, M ^ms + mc + mb. (13) 

Concurrently, the six-particle phase space of the final state of process (1) reduces to the 
four-particle phase space of the process 

e~{Pe-) + e+{Pe+) ^scbip) + s{Pi) + c(P5) + KPe), (14) 

while the probability of bound-state formation is controlled by the value of the baryon wave 
function at the origin of coordinates, the only model parameter in this approach. Eventually, 
the differential cross section for process (14) assumes the form 



(27r)4|M|2 |^(0)|%4, 

d^PA d^p5 d^pe d^p 

^ (27r)32E4 ' (27r)32E5 ' (27r)32£6 ' (27r)32£;' ^ ^ 

In evaluating the cross sections in question, we employed codes for integration that enter 
as ingredients into the CompHEP package [15]. As a necessary test, we first of all made sure 
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that numerical values of the cross sections are identical for different choices of the 4- vectors 
Os(T, Occr, a-bcT , and Oeo- , which are involved in the construction of the basic amplitudes. Having 
proven this, we prescribed ten iterations for the cross sections, each involving 100 000 steps 
of a Monte Carlo sampling of the integrand. The error in evaluating the total cross sections 
was 2.0-2.5 %, while the error in the differential cross sections was about 10 %, on average. 

Among theoretical uncertainties that affect cross-section values, the choice of renormal- 
ization scale in the running coupling constant for strong interaction, the values of the baryon 
wave function at the origin, and numerical values of the quark masses are of greatest im- 
portance. In what is concerned with the quark masses, the results of the calculations are 
the most sensitive to the choice of value for the lightest quark (strange one), because, for 
some gluon propagators, the minimal values of the denominators are 4m^. To illustrate this 
dependence, we everywhere present the results obtained at two values of the strange-quark 
mass, nis — 300 and 500 MeV. The remaining parameters were set to the following values: 
rric = 1500 MeV, nib = 4800 MeV, a = a{Mz) = 1/128.0, a., = a,(Mz/2) = 0.134, and 
sin^ = sin^ 9yf{Mz) ~ 0.2240; the numerical value of the wave function for the spin-3/2 
Qsc& baryon at zero relative coordinates of its quarks was borrowed from [16]: 

|V'(0)|2 = 0.90-10-=^ GeV^ (16) 

We note that the change in the characteristic energy scale as{^) in the running coupling 
constant Q;s(/i) from ^ = Mz/2 to fi = Mz leads to a change in the calculated cross sections 
by a common factor of [as{Mz)/as{M'z/2)]'^ = 0.665. 

For electron-positron collisions at ^/s = 91.2 GcV the table presents the values of the 
total cross sections ctot and the forward-backward asymmetry at the Z-boson pole. This 
asymmetry is defined as 

AFB^{crF-crB)/{crF + crB), (17) 

where apicTB) is the cross section for the production of ilgcb baryons traveling in the forward 
(backward) direction with respect to the direction of the electron momentum. 

Table 1. Features of flscb baryon production in electron-positron collisions at the Z-boson pole 



TTis, MeV 


o-tot, fb 


Afb 


300 


0.0534±0.0014 


0.162±0.024 


500 


0.0153±0.0004 


0.158±0.016 



Figure 2 displays the transverse-momentum (pt) and rapidity (Y) distributions of ^scb 
baryons at the strange-quark-mass values of — 300 and 500 MeV. For both values of 
the mass m^, the differential cross sections da/dpr peak at pr values approximately equal 
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to one-fourth of the total energy of colhding particles, while the quantities da/dY peak at 
small positive values of the rapidity Y. 

By using the concept of a fragmentation function, we can represent our numerical results 
in a simpler analytic form that is convenient for phenomenological applications. It is natural 
to break down the entire set of diagrams considered here into three groups that correspond 
to the fragmentation of b, c, and s quarks (in accordance with the flavor of quarks that 
are produced at the 7/Z vertex). Here, the fragmentation of b quarks plays a dominant 
role, whence it follows that, to a high precision, we can approximate the differential cross 
section as (see, for example, [4] and the discussion on the treatment of experimental data on 
electron-positron annihilation in [17, 18]) 

da/dz^a^b- Db-.n^^,{z), (18) 

where (j^j, is the total cross section for the process e~e~^ bb, Db^n^^S^) is the function that 
describes the fragmentation of a 6 quark into an Vtscb baryon, and the variable z is expressed 
in terms of the energy E of the final hadron and its longitudinal momentum p\\ as 

^ = (£; + P|,)/(^+P||)^ax. (19) 

For reasons of practical convenience, the variable z is often replaced by the variable 
— p/Pmax [17-20], which is close to it, or by xe — E/E^^ [21]. The distinction between 

these definitions vanishes in the limit of ultrahigh energies, but it can be sizable under actual 

conditions. 

Experimental results obtained for electron-positron annihilation are usually contrasted 
against the Peterson fragmentation function [22] 

~ .[1 -(!/.)-. /(I -.)P - 

where £ is a phenomenological parameter. 

If one disregards the aforementioned small asymmetry in the angular distribution of Vtgcb 
baryons and sets z xe, the relation between the differential distribution of the cross section 
with respect to the transverse momentum of the product baryon and the fragmentation 
function assumes the form 

da ^abbPT } Db^n^Jz)dz ^^^^ 



But if the variable Xp is used instead of it is necessary to set M = in relation (21). 
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The conclusions drawn from a comparison of relation (21) with our numerical results 
are as follows: if z ~ x^; , the parameter values of e = 0.098 ± 0.012 and 0.132 ± 0.018 
correspond to the strange-quark masses of = 300 and 500 MeV, respectively; if z ~ , 
the corresponding parameter values are e = 0.108 ± 0.016 and 0.147 ± 0.022. 

For the sake of comparison, we present values of the parameter e in the Peterson fragmen- 
tation function (20) that were obtained in experiments where electron-positron annihilation 
was explored at y/s = 10 GeV: e = 0.23to;Q6 for c-quark fragmentation into [17]; e = 
0.29 ± 0.06 for c-quark fragmentation into [19]; and e = 0.24 ±0.08 and e = 0.23lo:o6 for 
c-quark fragmentation into Sc according to the results obtained in [18] and [20] , respectively. 
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Fig. 2. Differential distributions of the cross section for ilscb baryon production in electron-positron 
collisions at the Z-boson pole with respect to the transverse momentum px (top left, bottom left) and the 
rapidity Y (top right, bottom right) at the strange-quark-mass values of = 300 MeV (top left, top right), 
and m-s = 500 MeV (bottom left, bottom right). Points represent the results of our Monte-Carlo calculations. 
The solid curves correspond to the calculations by formula (21) with the Peterson fragmentation function 
whose parameter takes the value of £ = 0.098 for = 300 MeV and the value of £ = 0.132 for = 500 
MeV. 
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